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TECHNICAL MEMORANDUM X- 73354 


LASER DOPPLER SYSTEMS IN ATMOSPHERIC TURBULENCE 


I. INTRODUCTION 

Marshall Space Flight Center (MSFC) has been developing laser Doppler 
systems, pulsed and continuous wave, for remote measurement of wind speed 
and turbulence for over a decade. A discussion of these systems is given by 
Huffaker [1-4] . Tte pulsed system, designed to operate over long atmospheric 
paths, is affected by attenuation due to atmospheric gases, particulate matter, 
and turbulence. The transmission loss experienced by the pulsed system has 
been worked out in detail in Reference 5. The effects of atmospheric turbulence 
on the performance of the laser Doppler system are considered in this report. 

A discussion of the basic effects experienced by a laser beam propagating in 
the normal atmosphere is presented first. 

The atmospheric turbulence manifested in random temperature fluctua- 
tions causes optical effects long Imown to astronomers as twinkling (variation 
of image brightness) , quivering (displacement of image from normal position) , 
tremor disk (smearing of the diffraction image), dancing (continuous move- 
ment of a star image about a mean point) , wandering ( slow oscillatory motions 
of the image with a period of approximately 1 min and angular excursions of a 
few seconds of arc) , pulsation or breathing (fairly rapid change of size of the 
image), image distortion, and boiling (time-varying nonuniform illumination 
in a larger spot image) [6] . In the case of a laser beam, the turbulence 
effects are broadly categorized as beam wander, spot dancing, beam spread, 
and scintillations. These effects are observed due to the temporal and spatial 
fluctuations of the direction, phase, and intensity of the wavefront. A descrip- 
tion of the atmospheric inhomogeneities causing the effects mentioned is given 
in the next section. 


1 1. PARAMETERS OF ATMOSPHERIC TURBULENCE 


The optical effects of interest are produced by the variations of refrac- 
tive index along the path of the beam. The changes in the refractive index are 
caused by the fluctuations of temperature which arise in turbulent mixing of 
various thermal layers. It is found that the temperature fluctuations obey the 
same spectral law as the velocity fluctuations. In analogy with the velocity 
turbulence, the atmosphere may be imagined to consist of a large number of 
eddies with varying dimensions and refractive indicies. For isotropic tur- 
bulence, the velocity spectrum (where the wavenumber K = 27r/f , i 

being the size of a turbulent eddy) contains the information about the turbulence. 
0 ( K) is intuitively interpreted as the amount of energy in the turbulent eddies 

of size i . The kinetic energy of turbulence is assumed to be introduced through 

eddy scale sizes larger than the ’’outer scale” of turbulence, Lq, corresponding 

to a spatial wavenumber Kq = 27t/ Lo. As the wavenumber increases beyond Kq, 

the turbulence tends to become isotropic and homogeneous. Experiments have 

confirmed the predictions of Kolmogorov* s theory of dependence for the 

three-dimensional spectrum up to wavenumber K related to tlie inner scale 

m 

of turbulence i o = 5. 92/f q* The spectrum in the dissipation region for 

which the scale sizes are smaller than n or wavenumbers greater than K is 

m 

steeper than The spatial wavenumber region between the inner scale 

£ 0 and the outer scale Lq obeying the Kolmogorov theory of a spectral slope of 
-11/ 3 is known as the inertial subrange, i q is of the order of several milli- 
meters and Lq several meters. 

At optical wavelengths the refractive index may be related to the tem- 
perature through the relation 



where n is the refractive index, P is pressure in mb, T is temperature in K, 
and X is wavelength in pm. The refraetive index structure constant a 

measure of the strength of atmospheric turbulence and is defined as the 


2 


mean-square difference in the refractive index at two points divided by the 
separation distance r raised to the 2/ 3 power, i. e. , 


^2 ^ ni)^_> 

n j,^/3 


( 2 ) 


The temperature structure constant is similarly defined for tem- 
perature fluctuations. Thus may be determined from temperature measure- 
ments alone from the following relation; 


C 

n 


77.6 P 


. 0.00753 


] 


X 10”^ a 


( 3 ) 


Typical temperature variations of approximately 1 K resulting in refrac- 
tive index variations of a few parts in a million will cause significant effects on 
optical radiation fields propagating through the atmosphere. is measured in 

\mits of (meters) and varies from 10”^® or more for strong turbulence to 
10”*^ or less for weak turbulence. 

Assuming that the form of the temperature spectrum and the refractive 
index spectrum is the same as that of the velocity spectrum, Tatarski [7] gave 
the following form for the refractive index spectrum: 


(A (K) = 0.033 exp 

n n 



(4) 


where 0 ~ spectrum has a singularity at K = 0. Physically the 

singularity implies that the energy per unit volume becomes unbounded as the 
eddy size increases. This singularity is avoided by using the modified von 
Karman spectrum [8]; 
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which gives a bounded variation for K < 27 t/ Lq. The three dimensional refrac- 
tive index spectrum models according to equations (4) and (5) are shown in 
Figure 1. Deviations from these models occur if the turbulence is not homo- 
geneous. Deviations have been observed in measurements over paths close to 
the ground and under conditions of weak turbulence [8] . The optical effects 
of the eddies are considered in the next section. 


11 1. OPTICAL PROPERTIES OF THE TURBULENT EDDIES 


The optics of the turbulent atmosphere may be described in terms of a 
collection of weak, moving, three-dimensional, gaseous lenses of varying scale 
sizes bounded by the inner scale i o the outer scale Lo of the turbulence. 
From equation (2) , it may be seen that the refractive index fluctuations asso- 
ciated with any scale size £ increases as the size of the eddy increases; 

<[6n(f)]^>= . (6) 


An eddy of radius £ can act as a lens with a focal length f = R/n^, where 
R is the radius of lens curvature. This is approximated as 

f ~ ~ 

n| n 

from equation (6) . Turbulence as it naturally occurs can consist of eddies of 
hot and cold air. Eddies cooler than the ambient air have a higher index of 
refraction. Li this case, nj is positive and the lens is convergent. K the eddy 
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SPECTRUM AMPLITUDE 



K - 27T/£ where £ - SIZE OF TURBULENCE EDDY; 
Kq - 2jt/Lq, typically Lq - 1 - 100 m; 

- 27T/Lfl, TYPICALLY £q - 1 - 3 mm 


Figure 1. Three-dimensional turbulence spectrum models 
for index -of-refraction fluctuation. 
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consists of warmer air, the index of refraction is lower (giving negative nj) 
and the lens is divergent. Typical values in the atmosphere are ~ 10”^, 

|nil ~ 10 “®, and i a few centimeters. Thus the focal length of the atmospheric 
lenses isf ~ 10 km. 

The turbulent eddies cause diffraction and refraction of the optical 

wavefront, and it is instructive to consider these two effects separately. Each 

eddy of size f bends a ray bundle individually and has a diffraction pattern of 

approximate width, \L/i , where L is the distance from the eddy to the receiver. 

The diffraction effects are negligible if XL/ S. « Ji or if i » n/ \L,' and this is 

the geometric optics regime. But if XL/ S. » i the diffraction pattern entirely 

determines the intensity distribution of the receiver. Thus the distance L , 

d 

beyond which the diffraction effects are important may be defined by ~ i^/ X. 

The refractive bending of rays by an eddy of size Jt depends on the 
refractive index of the eddy and may be calculated by Snell* s law to give 6 6 ~ 

6n. 6n is zero on the average. The mean square angular fluctuation for one 
eddy is <(60)^> ~ <(6n)S. 


Over a path length L, let the beam traverse L/ i eddies. Since the 
deflections are uncorrelated from one eddy to another, they add up as mean 
squares and the total mean square deflection is 


<<>~7 


( 7 ) 


The effect of refraction is to change the cross section of the beam from 
to {i + LO )^. The amplitudes of the wave before and after refraction are 
related as 


= (A+AA)2 (i ± 


The fractional amplitude change is then given by 


AA 

A 


L 


e 

r 


X 


/sj 



The variance of the amplitude fluctuation is then given by 


<x^> ~ F <«;> ~ 


( 8 ) 


The particular scale size 1 which contributes predominantly to the amplitude 
fluctuation depends on the path length; i. e. , different eddies are more effective 
at different path lengths in producing amplitude fluctuations. This point has 
been discussed by Tatarski [9] and De Wolf [10] as follows. 

There are no diffraction effects in equation (8) . Diffraction effects are 

negligible if the path length L is short or L « the geometric optics 

s s 

regime. Since I q (biner scale length) is the smallest eddy size, the diffraction 
by eddies of all sizes is negligible in this case. The smallest eddies then have 
the largest contribution. Equation (8), then becomes 

<X^> ~ for . (9) 


For path lengths greater than L , such that L » il/X, the eddies 
£ S i 

with sizes greater thani^ (f > have negligible diffraction, while eddies with 
sizes smaller than (i < diffract the rays where 

= n(xl7 

^ Jl 


The diffraction pattern entirely determines the intensity distribution and the 
refraction is negligible for the small eddies of size £ « ; i.e. , these eddies 

simply scatter the energy, not focusing the rays, and have little effect on the 
amplitude fluctuations. Geometric optics results may still be applied if the 
eddies of size £ < £^ are excluded. Then the smallest eddy size contributing 

significantly to the amplitude fluctuation is and ^ is substituted in 
equation (8): 
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( 10 ) 



n i 


-Ve 


n i 


for L, » 


b. 


This is also the result obtained by the Rytov approximation. Equation (lo) has 
been interpreted by Tatarski [9] as the result of the application of geometrical 
optics to a medium from which all eddies of scale sizes less than ^/\l7have been 
removed. 

A qualitative description of the effect of the optical inhomogeneities on the 
wavefront is sometimes helpful and is presented here. As a plane wavefront 
passes through a region of relatively lower refractive index, the wavefront 
will advance as the speed of light is higher. When passing through a region of 
higher refractive index, the wavefront is slowed down. Thus the refractive 
inhomogeneities of the atmosphere deform the wavefronts. Portions of the 
wavefront may be convex or concave in the direction of propagation. The con- 
cave portions converge while the convex portions diverge. Thus it is possible 
for the intensity to build up in some areas. Since the inhomogeneities are 
random and moving, the bright and the dark areas move randomly. This point 
of view has been developed by the optical astronomers [6] . 

Young [1.1] obtained equation (lo) on considering an aperture of 
Fresnel-zone size while calculating the normalized variance or the total modula- 
tion power obtained with a circular aperture using geometric optics. He 
observed that ray optics gave the correct parametric dependence because very 
little power is contributed by the range of spatial frequencies greater than 
-(\L) which are in the wave optics regime. Thus weak scintillation may be 
thought of as aperture filtering with a Fresnel-zone size aperture which filters 
out spatial frequencies higher than 

A related view developed by the radio astronomers is the phase screen 
concept. Immediately after passing through a layer of inhomogeneities, only 
the phase is affected and the intensity remains constant. As the distorted 
wavefront propagates further, the relative phases will be changed due to 
different directions of propagation and in this way, the amplitude fluctuations 
develop. This concept has been applied by Lee and Harp [12] to wave propaga- 
tion by dividing the three-dimensional refractive field of the medium into thin 
slabs perpendicular to the direction of propagation. The effects of the refrac- 
tive inhomogeneities on the propagation of laser beams in the atmosphere is 
discussed in the following section. 
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IV. EFFECT OF ATMOSPHERIC TURBULENCE ON LASER BEAMS 


A qualitative description of the effects based on the relative sizes of 
the beam and the refractive inhomogeneity is given first, followed by the results 
of Rytov approximation on the wave equation. 

If the scale size of the inhomogeneity is much larger than the diameter 
of the laser beam, the entire beam is bent away from the line of sight and 
results in beam wander or beam steering. The center of the beam executes a 
two-dimensional random walk in the receiver plane and, for isotropic turbulence, 
the displacement of the beam from the line of sight will be Raleigh-distributed [ 131 • 
Inhomogeneities of the size of the beam diameter act as weak lenses on the whole 
or parts of the beam with a small amount of steering and spreading. When the 
inhomogeneities are much smaller than the beam diameter, small portions of the 
beam are independently diffracted and refracted and the phase front becomes 
corrugated. The propagation of the distorted phase front causes constructive 
interference over some parts of the receiver and destructive interference over 
others, leading to alternate bright and dark areas. Since the atmosphere is 
seldom stationary, the locations of the bright and dark regions change con- 
tinually and give the pattern of ”boiling." Thus, the turbulence causes the 
received field to scintillate in time and space. Since the atmosphere consists 
of inhomogeneities of all sizes in motion, the laser beam experiences fluctua- 
tions of beam size ( spreading) , beam position ( steering or wander) , and 
intensity distribution within the beam ( scintillation) simultaneously. The 
relative importance of these effects depends on the path length, strength of 
turbulence, and the wavelength of the laser radiation. The variance of irradi- 
ance reflects the deviations from the mean due to beam wander, broadening, 
and scintillations. Estimates of these effects for finite beams are discussed 
as follows. 

Finite transmitter effects have been treated by the Huygens -Fresnel 
approach [14] and by the transport approximation, Fante [15]. The Huygens- 
Fresnel principle was extended by Lutomirsld and Yura [14] so that the second- 
ary wavefront is determined by the envelope of spherical wavelets from the 
primary wavefront as in vacuum, but each wavelet is determined by the propaga- 
tion of a spherical wave in the refractive medium for small scattering angles. 

Hence a knowledge of the spherical wave propagation in the turbulent medium 
is sufficient to determine the field distribution at an observation point. The 
mean irradiance distribution at the observation point is characterized by the 
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mutual coherence function (MCF) of a spherical wave and the geometry of the 
transmitting aperture. The MCF, defined as the cross-correlation function 
of the complex fields in a direction transverse to the direction of propagation, 
describes the loss of coherence of an initially coherent wave propagating in a 
random medium. When the medium is characterized by a refractive index 
variation that is gaussian with zero mean, the MCF is given by [7] 


M(p,L) = exp 



( 11 ) 


where D is the wave structure function, p is the lateral separation between the 
two points in a plane perpendicular to the direction of propagation, and L is the 
path length of propagation. For spherical wave propagation and a Kolmogorov 
spectrum of turbulence [16] , the wave structure function is given by 


D(p) = 2.91 



For constant C^ , 
n 


D(p) = 1.089 K^C^Lp 



( 12 ) 


(13) 


Equation (ll) may be written as 


M(p, L) = exp 


[■( 5’1 ■ 


io « p « L, 


0 » 


(14) 


where 


Po = 


1.455 k 2 / dstf (s) 


I ”3/5 
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For constant 

n 


Po 



[0.545 K^LC^] 
n 


( 15 ) 


Po is the lateral separation such that the MCF becomes equal to 1/e and is 
called the lateral coherence length. If d > po where d is the diameter of the 
aperture, the turbulence along the path of the beam reduces the lateral coherence 
between different elements of the aperture and effectively transforms it into a 
partially coherent radiator of dimension, po, which decreases with increasing 
distance from the aperture [13] . But if d < po, the entire aperture behaves like 
a coherent radiator. Thus the resulting field distribution is cliaracteristic of a 
coherent aperture with a diameter d equal to tiie smaller of po and d. Yura [ 17] 

showed that if L » Kd th© beam properties may be approximated by those 

of a spherical wave and if L « Kd the beam properties may be approxi- 
mated by those of a plane wave. 

In the absence of turbulence, a gaussian laser beam will have a far-field 
angular spread 0 = 2\/7id. When turbulence is present along the path, the 
scattering of the beam by the moving turbulent eddies causes additional spread- 
ing which is greater than 0 under some conditions. 

The spreading of the beam consists of the deflection of the beam as a 
whole by the eddies larger than the beam diameter d and by the broadening of 
the beam due eddies smaller than d. The separation of the beam spread into 
the two components of wander and broadening is usually based on the length of 
the time of observation. If V is the transverse wind-velocity component, the 
turbulence eddies are interacting with the beam in time intervals of order 
At ~ d/V. For times much longer than At, the received spot will be fully 
spread due to wander and broadening. The total spot size, \dewed on a time 
scale T » At, is called the long-term spot size and the beam radius of long- 
term Spread is denoted by p If a picture with an exposure time much loss 

Xj 

than At were taken of the beam, a laser spot of radius p broadened by liie small 

s 

eddies at some distance p normal to the line of sight would be observed, p is 

C o 

the radius of beam-centroidal motion. 
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The concepts of the beam spread just discussed apply only if the beam 
is intact. When the beam wander effect is removed, the atmospheric spreading 
of the beam at the receiver is the short-term beam spread; however, when the 
turbulence along the path is strong, the beam breaks up into multiple patches 
and there is not much wander. Kerr and Dunphy [18] found experimentally at 
near field that the focused beam broke up at the receiver into a proliferation 
of transmitter-diffraction-scale patches when there is either transmitter mis- 
adjustment or strong turbulence. Reidt and Hohn [19] observed in recent 
experiments that the focused beam at 0. 63 pm was broken into several spots of 
characteristic length 2\L/7rd which is the free-space diffraction-limited focused 
radius. They verified also the inverse variation of the patch size with the 
diameter of the beam, that is, the increase in patch size with decreasing beam 
diameter or vice versa. 

The beam will be quivering while remaining intact as long as the tur- 
bulence is relatively weak and the path is not too long or, more precisely, if 

L < L where L is the critical range defined by K d^ where d 

CR CR CR eii Gii 

is the smaller of po and d [20] . 

Quantitatively, the spreading of the beam depends on the lateral coherence 
length pq. If the beam diameter d is much smaller than po, the effect of turbu- 
lence on the beam is negligible and the beam diameter in the far field will be 
that due to diffraction, namely, 2XL/7rd. However, if the beam diameter is 
greater than the lateral coherence length po, the beam diameter is given by 
2\L/7rpo [20] . That is, the beam behavior is characteristic of an aperture of 
diameter po instead of d. The turbulence along the path transforms the aperture 
into a partially coherent radiator. For the initial field distribution of gaussian 
form given by 



where F is the radius of the curvature, the long term beam radius is obtained 
as [21] : 
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When the turbulence is strong, the beam will be broken up and for this case 
Klyatskin and Kon [22] found that the wander of the center of gravity of the 
beam may be obtained from 

«> “ (17) 


for constant C^. L » Kjpo (far-field case) , and d » pQ. Taking the ratio of 

equations (17) and (16) and observing that the last term on the left hand side 
of equation (16) is dominant for strong turbulence, we obtain 



Since L » Kpo then \ L/ * means that the motion of the beam 

centroid is negligible compared with the beam spread when the beam breaks 
up into multiple patches. It is not possible to tell how many patches will form, 
but the bright patches will be in a zone with mean square radius [23]. 


For ranges less than the critical range, tracking laser transmitters have 
been proposed and discussed [24,25] to compensate the beam wander which is a 
relatively low-frequency phenomenon of the order of 0. 5 Hz. For L < the 

beam centroidal motion is given approximately by [26] 



2. 97 1 } 


(19) 


and appears to satisfy a gaussian distribution [27] on each axis. 
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One of the basic observations on laser beams is the fluctuations in 
intensity in addition to fluctuations in beam size and position. Fluctuations 
in intensity are called scintillations. Experimental data on scintillations are 
extensive at present and are usually expressed as the normalized variance of 
the intensity: 


w - (I>^ 

(l) ^ 


( 20 ) 


I is the measured value of the fluctuating intensity and (0 is the average. 
The statistical theories of optical propagation in a randomly inhomogeneous 
medium usually calculate the variance in the log-amplitude X. The log- 
amplitude X(x) at any location x is related to the Intensity or the irradiance 
I(x) by 


x(x) = In 1 In ^ (21) 

where A is the amplitude. The variance of the log-amplitude is theoretically 
given by 


0-2 - 0. 124 

T n 


L » 



( 22 ) 


for a spherical wave. Experiments have confirmed the increase of the log- 
amplitude variance as the 11/ 6 power of the path length, as the 7/ 6 power of 
the optical wavenumber K, and the square of the strength of turbulence for a 
homogeneous path, i. e. , = constant along the path as long as <0.3 for 

the spherical wave [28]. 

Increase of the path length or the strength of turbulence beyond cr^ = 

0. 3 fails to produce corresponding increase in the measured value of the log- 
amplitude variance cr^ which slows down, reaches a maximum, and begins 
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to decrease, reaching an as 5 nnptotic value of approximately 0.23 independent 
of wavelength, path length, and parameters of the turbulent medium. Thus the 
scintillation of an optical beam does not increase indefinitely but tends to 
saturate for strong turbulence and long paths. This is shown schematically in 
Figure 2. The more precise way of showing the saturation effect is to plot the 
theoretical value of log-amplitude variance along the x-axis and the experimental 
value of the log-amplitude variance on the y-axis, as is usually done in the 
literature [28]. 


The probability distribution of the amplitude or intensity fluctuation is an 
important aspect in optical communications. Theoretically the application of 
the central limit theorem leads to a normal distribution of the log-amplitude; 
i.e. , the scintillations follow a log-normal distribution. For a unit-amplitude 
plane wave, the probability density of intensity p(I) satisfies 


p(l) 


exp 

(27t) 


In 1 + 



where 


= In + 0 - j j 

For « 4, there is extensive experimental evidence in support of log-normal 

distribution of intensity. Physically the reason for the log-normal distribution 
is the multiplicative effect of the refractive -index variations along the path 
on the intensity [9] ; i. e. , the field is modulated in each subrange of the path 
multiplicatively if it is assumed that the effect of atmosphere in each subrange 
is independent of the degree of coherence. Hence the variation of the logarithm 
of the amplitude and intensity is the sum of several random variations induced 
along the path, of propagation which leads to a normal distribution due to the 
central-limit theorem. Experimental measurements by Gracheva, et al. [29] 
indicate that the log-normal distribution is a reasonable approximation for 
25 < o-^ < 100. For 1 ^ <7^ 25, there are significant deviations from the 
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LOG-AMPLITUDE VARIANCE 



Figure 2. Saturation of the variance of log-amplitude fluctuations of a spherical wave 
propagating through a homogeneous (C^ = const.) turbulence path. 


log-nonnal distribution. They did not notice Rayleigh distribution, as predicted 
by de Wolf [30] for large values of Physically, Rayleigh distribution 

results when the field at the receiver is the sum of a large number of independ- 
ently scattered fields. The physical model of de Wolf [ 10] attributes the 
Rayleigh distribution to the multiple scattering by off-axis eddies which are 
large in strong turbulence. 

In addition to temporal variation, the received field exhibits a spatially 
varying intensity after propagation through turbulence. The spatial structure 
of scintillations is usually studied by the covariance of the logarithm of the 
amplitude fluctuations. A correlation length may be defined in several ways. 

One definition is the value of separation for which the spatial covariance ftmction 
goes to zero for the first time. The correlation length provides a measure of 
the separation between two points over which the scintillations are correlated. 
The concept of the correlation length has been useful for some applications. 

For example, the theoretical calculation of the log-amplitude variance is 
usually based on a point detector. In practice, a collector aperture of diameter 
less than the correlation length is treated as a point detector for comparison of 
theory and experiment. Another use of the correlation length is in understanding 
the aperture averaging effect which is discussed later. The intensity covariance 
function is usually measured and is related to the log-amplitude covariance 
as [31] 

Cj(p) = <l) [exp{4C^(p)} -l] (23) 

where the separation is p = Ix - x* |. the covariance of log-amplitude is 
C^Cp) = <[X(x) - (x>) [X(x') - (x)l> 
and the covariance of intensity is 

Cj(p) = (ir(x) - (i)j [I(x>) - (I>]> 
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The intensity and the log-amplitude correlation distances may be seen 
to be the same from equation (23) but the curves differ greatly. For plane or 
spherical wave propagation, C^(p) the form 


C^(p) 



where f is a function which depends on the mode of propagation. The first zero 
of f is the correlation length and is approximately 0. 76 for a plane wave 

and 1. 8 (\L) for a spherical wave based on the perturbation theory of 
Tatarski. Experiments confirm a correlation length of the order of (\L) 
when there is no saturation of scintillations. 

The classical method of reducing the fluctuations in the received power 
has been to increase the size of the receiving optics. The reduction of the fluc- 
tuations in the received power, when the diameter of the receiving aperture is 
Increased, is called aperture averaging. The physical reasoning is that large 
apertures tend to average over the statistically independent portions of the 
scintillation pattern [ 31] , The diameter of an independent patch is usually 
imagined to be of the order of the correlation length. Thus for weak scintilla- 
tions, the size of the patch is of the order of the Fresnel length Experi- 

ments [32] show that the correlation length decreases as the turbulence gets 
stronger. Moreover, the covariance curve develops a progressively higher 
correlation tail and aperture averaging deteriorates as the turbulence strength 
increases. Figure 3 illustrates the covariance function for several o-^. 

To summarize, the phenomenon of saturation of optical scintillation 

■ . t 

may be characterized by the following experimental observations: 

1. There is a physical break-up of the beam into multiple patches when 
the turbulence is strong. 

2. The variance of the log-amplitude reaches a maximum instead of 
increasing indefinitely and then decreases to an asymptotic value of approxi- 
mately 0.23 independent of the wavelength of the optical radiation and param- 
eters of the turbulent medium as the path length, or the strength of turbulence, 
or both are increased. 
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Figure 3. Theoretical curves of the log-amplitude covariance function in weak 

to strong turbulence [ 34] . 


CO 


3. The spatial covariance function of log-amplitude decays progressively 
faster but maintains a progressively higher tail as a function of the detector 
spacing with increasing turbulence strength and path length. 

4. Aperture averaging becomes progressively less effective for a given 
size of the aperture as the turbulence strength and the path length are Increased. 

There were several theoretical attempts to explain the saturation 
phenomenon in the past but it is only in recent years that the theory is able to 
accoimt for all the observed effects. Since the saturation effect is observed 
under conditions of strong turbulence, a working definition is needed for 
the conditions specified as weak or strong turbulence. The turbulence is 
called "weak" if the theoretical variance of log-amplitude 

0-^ = 0.124 
T n 


satisfies the inequality cr^ < 0. 3. The turbulence is called ’’strong” if cr^ > 0. 3. 

Thus the specification of the turbulence conditions involves not just the refractive- 

index structure constant but a combination of the wavelength of radiation, path 

n 

length and C^. 

One of the earlier attempts to physically explain saturation is that by 
Young [33] . When stellar scintillation is observed with small apertures, the 
total modulation power first increases with zenith angle and then decreases. 

Yoimg suggested that the smearing or blurring of the details of the wavefront 
distortions caused by the intervening atmospheric turbulence is responsible 
for the observed saturation of astronomical scintillation. The blurring effect 
is produced by the small-scale fluctuations in the wavefront which change 
rapidly due to the short lifetime of the small eddies. Thus, a shadow pattern 
produced by a high layer in the atmosphere becomes smeared or washed out by 
the intervening turbulence as it propagates down through the atmosphere. 

Clifford, Ochs, and Lawrence [34] give a detailed description of the 
smearing process in the scintillation pattern in strong turbulence. The wave- 
front develops small-scale and large-scale distortions while passing through the 
turbulent medium. When such a distorted wave is incident on a Fresnel-zone- 
size eddy, the small-scale distortions reduce the resolving power of the eddy 
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while distortions larger than the eddy merely tilt the diffraction pattern of the 
eddy. The small scale distortions may be expected to change in detail several 
times during the lifetime of the larger eddies due to the shorter lifetimes of 
small eddies. These fluctuations of the small-scale distortions of the wavefront 
cause the smearing of the diffraction pattern of the Fresnel-size eddy, thus 
making it less effective as a producer of scintillation. Scintillations remain 
bounded and do not build up indefinitely due to the smearing by the small-scale 
eddies. 

The smearing effect of strong turbulence is Incorporated into the first- 
order analysis by means of a spectral -filter function which filters out higher 
spatial frequencies. The spectral filter function is given by the normalized 
two-dimensional Fourier transform of the short term average of the irradiance 
profile of an optical beam propagating through turbulence. The filter function 
turns out to be the short term modulation transfer function (MTF) . It is 
assumed that the log-amplitude covariance of a spherical wave for strong 
turbulence is given by the first-order log-amplitude covariance whose spatial 
frequency spectrum has been modified by multiplication.with the spectral filter 
function (or the short term MTF). The log-amplitude covariance thus calculated 
agrees well with the experiments and is shown in Figure 3. 

From this covariance function, Clifford and Yura [35] obtained the 
follo’w'ing asymptotic behavior for the log-amplitude variance: 

al. = 0.36 (o-2) ^ 0.42 (24) 

JL 


where o; is a constant found from experiments. For a = 0.7, which is obtained 
from the strong scintillation data, ar^ tends to 0. 23 for large values of o-^. 

The perturbation theory of Tatarsld [7,9] implicitly assumes that the 
Initial coherence is maintained along the propagation path and eddies having 
scale sizes of the order of the Fresnel length are most effective in producing 
scintillation. However, as the wave travels through the medium, the transverse 
spatial coherence decreases continuously and the wave becomes partially 
coherent. As was previously noted, the spreading of a beam of initial diameter 
d in a turbident medium is characteristic of an aperture of diameter po (the 
lateral coherence length) and not of diameter d. The radiator has become 
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partially coherent in the presence of turbulence. Incorporating the loss of 
coherence of a radiator in the turbulent medium into the analysis, Yura [36J 
obtained the saturation of scintillations and the behavior of the amplitude 
covariance. He showed also that in the strong turbulence regime, the ampli- 
tude correlation length is equal to the lateral coherence length and increase. 

The turbulent eddies most effective in producing scintillations in the saturation 
regime are those that have scale sizes of the order of the lateral coherence 
length and not the Fresnel length. 


9 "Vs 

The asymptotic log-amplitude variance behavior of (o-^) has been 

predicted also by Fante [37] and Gochelashvily and Shishov [38] from approxi 
mate solutions to the equation satisfied by the fourth moment of the field. 


The physical theories [34, 35, 36] and the approximate solutions [37, 38] 
have been able to predict the observed behavior of the covariance function, 
namely, a rapid fall -off followed by a slow decay resulting in a long tail. These 
curves indicate that there are small-scale and large-scale structures in the 
scintillation pattern for strong turbulence. The correlation in strong turbulence 
is over a transverse separation of the order of 


Po 


(Mi 


V2 




V5 


As increases, the lateral coherence length and the amplitude correlation 

length decrease. The diameter of the detector must be smaller than po for 
point detector performance in strong turbulence. 

The aperture averaging which produced a reduction in scintillat" on noise 
became less effective for a given aperture of the receiver as the turbulence 
increased. This effect is due to the presence of large-scale scintillations at 
the receiver as o-^ increased. 

Another quantity of interest is the temporal frequency spectrum of 
intensity fluctuations observed at a fixed point. The width of the frequency 
spectrum is of the order of V/ (xL) for weak turbulence where V is the 
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velocity component normal to the beam. Tliis is typically 10 to 100 Hz. For 
the case of strong turbulence, the width of the frequency spectrum is of the 
order V/pg which is approximately 100 to 1000 Hz. In fact, the peak of the 
spectrum occurs at approximately V/(XL)*'^ for cr^ < 1 and at approximately 
V/po for cr^ > 1. 

In addition to the shifting of the peak of the spectrum to higher fre- 
quencies, there is a broadening of the spectrum and the maximum value of the 
spectrum decreases for Increasing cr^. These features of the spectrum have 

been qualitatively observed in the experiments. Finally, to obtain all of the 
temporal Information of the scintillation in the strong tui’bulence regime, it may 
be necessary to use an electronic bandwidth of the order of V/ pg. Otherwise, the 
use of a smaller electronic bandwidth may result in an apparent decrease of 
the log-amplitude variance for increasing values of or^. 

V. EFFECT OF TURBULENCE ON HETERODYNE 
LASER DOPPLER SYSTEMS 


In a heterodyne receiver, the incoming wave which is usually weak is 
combined with a reference wave on the photodetector surface. Coherent detec- 
tion results from the square-law response of the photodetector to the incident 
radiation. The squaring operation of the photodetectcr produces a current at 
the intermediate frequency which contains the signal -modulation information. 

The heterodyne detection of an atmospherically distorted wavefront 
has been the subject of several studies [39, 40] . The various effects of atmos- 
pheric turbulence on the laser beam are related to the image behavior in a 
heterodyne system. Comparison of the expression for the optical power density 
due to a laser transmitter at a point some distance L aw'ay from the transmitter 
and the expression for the signal-to-noise (S/N) ratio of a heterodyne receiver 
detecting a signal from a point source at a distance L shows that they have 
identical functional dependence on the turbulence parameters. Thus the tele- 
scope performance measured in terms of the coherence diameter remains the 
same whether it is used as part of a transmitter or as part of a heterodyne 
receiver for propagation in a turbulent atmosphere [41] . In simple terms, 
the reciprocity between the transmitter and receiver tells that the average target 


23 


illumination is related to the power received by the heterodyne receiver. Fur- 
ther, beam wander and beam spread in the target illumination system are 
related respectively to image dancing and spread in the heterodyne receiver. 
Thus the beam-image reciprocity implies the familiar effects of turbulence 
on the beam propagation in the analysis of the heterodyne detection. 

Fried [39] considered the amplitude and phase fluctuations and found that 
for an atmospherically distorted plane wave and a plane-wave local oscillator 
there is an upper limit to the achievable S/ N no matter how large the detector 
aperture is. The aperture diameter for the near -peak value of S/ N is of the 
order of the lateral coherence diameter. This aperture size has the virtue 
of avoiding large signal power fluctuations and variance which would result for 
larger collection apertures. Moreland and Collins [4uj tried to determine the 
possibility of increasing S/N by properly sliaping the local oscillator beam and 
came to the conclusion that, though the optimum local oscillator shape differs 
significantly from that of a plane wave, the increase of average S/N over the 
plane-wave case is negligible. 


VI. CALCULATION OF S/N LOSS 
FORTHEMSFC SYSTEM 


The MSFC pulsed laser Doppler system consists of a coherent CO 2 laser 
operating at 10. 6 fxm wavelength and transmits 140 pps adjustable in width from 
2 to 8 ixs through a modified Cassegrain telescope. A homodyne receiver collects 
the backscattered laser radiation which is detected by an infrared detector 1 
and processed in signal-processing electronics. A block diagram of the MSFC 
system is shown in Figure 4. The transmitted laser beam is as.sumed to have 
a gaussian amplitude distribution in the radial direction. The laser beam passes 
through a telescope with a nominal range of focusing f from the transmitter. 

The wave function at a distance L from the transmitting telescope of 
radius R may be written by the Fresnel diffraction theory in the absence of tur- 
bulence along the path as [42-44] 
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where 


!/'(x,y, L) 



[exp i(KL - cut) 1 
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II 
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dx" dy” 


( 26 ) 


= number of photons transmitted per second, 

r" = vector in the plane of the transmitting lens , 

r = vector in the plane at distance L from the transmitter and normal 
to the optic axis, 

A = wavelength of laser radiation, 

cu = angular frequency 

t = time. 

Equation (26) may be integrated by completing the square and the result 
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(27) 
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The effect of the atmospheric turbulence along the path of propagation is 
to produce random variations of phase and intensity. These random variations 
are introduced into the analysis by multipl 5 dng equation (26) with 


[X(r, r”) + i S(r, r”)j 


X( r , r") is the perturbation of logarithm of amplitude, equation (21) , and 
sCir , r”) is the perturbation of phase between the point at the transmitter 
and the range point r . Thus the outgoing wave function at a distance L with the 
effect of turbulence is given by 

/ 2V/2 A r% 

ip{x,y,L) = I -j exp [i(KL - cot)] f f exp 

' ’ —CO 

- i ^ + i S(T, T”) +X(7, T”) - -^1 dx” dy” . (28) 
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The transmitted wave is assumed to be backscattered by a diffuse target 
of scattering particles. Considering a single particle, the incoming wave is 
scattered by the particle as a spherical wave, and the scattered wave function 
at a point *r» in the plane of the receiving lens for a coaxial system is given by 



where 


(T^ = backscattering coefficient of the particle. 


A oj = 2u)V/ c 


Doppler shift due to the velocity component V of the 
moving particle along the optic axis. 


X(r, r') = perturbation of the log-amplitude, and 

S( r, r ' ) = perturbation of phase, between the range point r and the 
point ’r* at the receiver . 

The local oscillator reference wave function which is added to the 
received wave function is assumed a plane wave with gaussian amplitude varia- 
tion which may be written as 



( 30 ) 


where a is a constant and the arbitrary phase factor is assumed to be zero. 
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The total signal photo-current is given by the real part of 


1 = 2t] f if) ip* -W(r*) dr’ 

= ■’ s ref 


where W(r’) is the apodization or transmission function for the receiver 
aperture and is assumed equal to unity over the width of the reference beam. 

Substitution of equations (28), (29), and (30) into equation (3l) gives 
the signal current as the real part of 


ff exp -|•^ + •^| + i|Aa5t+(^)(r’^ r», r) 


• - 2(T^aP3 

^s "" ttlFrx 


+ S( r, r”)+S(r, r’ )/+ X( r, r”) + X( r, r’) dr* dr” 




where 


cp (7, 7* , 7”) = 2KL + ^ j(7 -7)2 + (7 -7”)2j- ~ (7*2 +7*2) 


The photo-current due to a single scatterer may be written as 


^s 7 tI?R\ 


2 j,tii2 1 ^ ^ 

r+-^( +X(r, r*) 


oiAhi r r f ( r*2 r*”2 \ ^ 

= ' // “"p [- 1-5^ j ^ 

+ X(r, cos ^Aa)t + ^(r", r*, r) + S(r, r*) 


+ S( r, r*») I dr* dr” 
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The signal power is proportional to i^ and 
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+ x(r, 7») + x(7, 7”) + X(7, 7»”) + x(7, r””)j cos j^i 
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■ ^ ////.» 

+ X(7, 7") + x(7, r»") + X( r , r"") 

+ 0(7",7,7) +<a( 7"",7",7) + s(7,7) + s(r, r") 

+ S(7, 7") + s(7, 7"")] + cos (7”, 7 , 7) + S( r , r») 

+ S(7, 7") -0(7""+7"+r) -S(r, r*") 

- s{7,7"") 11 dr* dr" dr*" dr"" 

If the signal power is averaged over a period longer than 1/ (2Aco) , then 
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where Re denotes the real part. The point r^” lies on the receiver aperture and 
"rtt” lies on the transmitting aperture. 


The expected value of the signal power is obtained by taking the ensemble 


average of i^ as given by 
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,2a4^2 
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+ i (0{ r”, r* , r) - r’’”, r*”, r)! 


exp 


X(r, r») 


+ X( r , r»”) + X( r , r") + X( r , r’"’)|+ i js( r , r*) + S( r , 7") 

+ S(7, 7»”) “ S(7, V dr» dT” d7»” dr*"’ . (33) 


It is assumed that the outgoing and incoming paths are statistically 
independent. This assumption allows the following: 



r*) + X( r , r’") / + i (s( r , r*) - S( r , r’”) 


+ {X(r, r”) +X(r, r"”)j +i|s(r, r”) 
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- S( r , r"”)) j > = exp 

- - 1 

+ D( I r" - r’"’|) 
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where D( r) is the wave structure function. Kolmogorov theorj'^ of isotropic 
turbulence gives a 5/ 3-power variation for the wave structure function: 
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where is proportional to the coherence radius and has the value 

6/5 “3/5 -6/5 

r = 0. 05972 \ L C , for plane wave 
a n 


= 0.1047 , for spherical wave 


It will be appropriate to use r^ for a spherical wave for the wave originating at 
the range point. 

To simplify the calculations, the structure function is approximated as 


■ (tj 


and the plane -wave value of r is used for both ways to give a conservative 
estimate. 

The expected value of i^ becomes 

s 


< 0 = 


cr* 


,//// 


exp 


j-f 2 ^ J,ff2 ^ ^ J,ttf|2 

r2 


_ +i U(r’», r», r) 


- <?i{r"", r”, r) 


dr» dr” dr»” dr’”’ 


(34) 


Performing the integrations approximately, the expected value of SNR 
for a single scatterer is obtained as [42] 
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exp 


(35) 




where a = (R/r ) and t is the pulse duration. Integrating equation (35) over a 
& 

plane normal to the axis at the range L, the S/N ratio due to a plane target 
(dropping the angle brackets) is obtained; 


S_ 

N 


2a2 


ncr^TR A 

l 2(1 + a^) + ^ (f - L)2 


(36) 


In the absence of the turbulence effects, a = 0 and thus the S/N is reduced by 
the ratio 


P = 


L^(l + a=') + ^(f-L)' 


(37) 


due to turbulence. If the system is focused to infinity, then 




p = 


L2(l+a2)+^ 

L* + 2fSl 

L + ^2 


(38) 


Thus the calculation of the loss of power due to turbulence reduces to 

mlation of a = R/ r and P. 

a 

The wave structure function for a plane wave is given by 


D(r) = 2.91 


f c2 dL 
^ n 
Path 


(39) 
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For a propagation path which is inclined at an angle 0* with respect to the 
horizontal, 


h = H + L sin 0’ 


( 40 ) 


where h is any altitude above sea level and H is the altitude of the scattering 
surface above sea level. For flights at Edwards Air Force Base, the dry lake 
is approximately 600 m above sea level. The strength of turbulence is assumed 
to vary as 


2 2 v-1/3 

C = C _(h*) exp 
n nO' ' 



(41) 


2 2 
where C „ is a constant, h» = h - H, and h„ = 3200 m. The value of C 

nO 0 nO 

obtained from earlier measurements of Huffnagel and Stanley [ 39, 45] , is 

4. 2 X 10”^^. We use H = 600 m and assume that equation (4l) is valid from 

h = 10 m. Substituting equation (41) into equation (39) , the structure function 

may be written as 

D(r) = 3.38xlo'‘c^jjr^^® (42) 

X 

where l(i ) is the incomplete gamma function given by 

r(|)i(i) =/u"^* "du . i = 

and and are dimensionless lower and upper altitudes. From equation (42) 
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r = 
a 


4 2 cosec 0* 

3. 38 X 10 C „ ■ 
nO 




-3/5 


For the flight tests, 0 = 7® and X = 10. G fim. Using these values, r 


becomes 


r = 5.8 X 10 
a 


-10 




-3/5 


(43) 


The optics diameter is taken as 25, 4 cm and this gives R = 0, 127 ni; therefore, 


a = ^ = 2. 19 X 10® I I(i ^) - I(i 
^a 




(44) 


The incomplete gamma functions may be obtained from the Pearson Tables [46] . 

The calculations are made for the flight tests bet\\'cen GOO m and 4000 m 
2 -14 

altitude and for * 4. 2 x 10 . The calculated S/N values shown in Figure 

5, also include the transmission loss taken from Reference 5. 

The loss due to the atmospheric turbulence is smaller than that due to 
absorption at 10. 6 nm wavelength. As observed in Reference 5, the absoiption 
loss will be least at the DF (deuterium fluoride) laser wavelength, 3. 9 jam, but 
the turbulence loss may be higher due to the smaller lateral coherence diameter. 

A comparison of the losses due to absorption and turbulence over a 20 1cm 
horizontal path at altitudes of 2. 5 Ion and 5 1cm is given in Table 1 at CO 2 and DF 
laser wavelengths using the same optic size. Due to the smaller coherence 
diameter at DF wavelength, the optics diameter should be less than that at CO 2 
wavelength. By using the same optics size for botli, the turbulence losses have 
been overestimated for the DF laser. Eyen so, the total loss for DF system is 
approximately 15 dB less than for CO" system at 2. 5 1cm. At '5 Icm altitude, the 
magnitude of the losses has dropped and there is an improvement of only 5 dB 
using the DF system. For altitudes below 5 1cm,. the DF system is seen to give 
less total loss due to the atmosphere. 
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Figure 5. Comparison of measured and theoretical S/ N values 



TABLE 1. COMPARISON OF LOSSES AT CO 2 AND DF LASER 
WAVELENGTHS OVER A 20 1cm HORIZONTAL PATH 
AT 2. 5 1cm AND 5 1cm ALTITUDES 


Altitude 

Laser 

Absorption 

Loss 

(dB) 

Turbulence 

Loss 

(dB) 

Total Loss 
(dB) 


CO 2 

24 

2.6 

26.6 

2. 5 km 

DF 

1.75 



9.05 

10.8 


CO 2 

9.6 

0.37 

10.0 

5 1cm 


0.46 



DF 

4.9 

5.4 


Since the atmospheric conditions are constantly changing, it is desirable 
to estimate the maximum and minimum attenuation of the CO 2 laser radiation 
that may be expected due to the atmospheiic propagation. The upper limit on 
the attenuation is estimated using the AFCIRL Midlatitude Smnmer Hazy Atmos- 
phere and an exponential variation with altitude of the turbidence strength witli 


2 -14 2/3 

C „ = 4. 2 X 10 m . The lower limit on the attenuation is estimated using 
nO 

AFCRL Midlatitude Winter Clear Atmosphere and an exponential variation with 


2 -15-2/3 

altitude of the turbulence strength with C^^^ = 4. 2 x lo m 


The results 


for a two-way horizontal path of 20 1cm range at each altitude up to 10 1cm 
altitude are shown in Figure 6. The difference in the upper and lower limits 
of attenuation is considerable near sea level mainly because of the differences 
in the temperature, humidity, and aerosol concentration. The brealaip of the 
attenuation into the transmission part and turbulence part at each altitude is 
given in Table 2. 
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ATMOSPHERIC ATTENUATION - dB 



Figure 6. Atmospheric attenuation at 10. 6 pm. 



TABLE 2. SIGNAL LOSS IN ATMOSPHERIC PROPAGATION 
AT 10. 6 fxm DUE TO ABSORPTION, SCATTERING 
AND COHERENCE DEGRADATION FOR A TWO-WAY 
HORIZONTAL PATH OF 20 km AT EACH ALTITUDE 


2 

(a) AFCRL Midlatitude Summer Hazy Atmosphere and C = 
4. 2 X 10-‘« 


Absorption and 

Turbulence 



Scattering Loss 

Loss 

Total Loss 

Altitude 

(dB) 

(dB) 

(dB) 

10,0 m 

68.0 

13.0 

81.0 

2. 5 km 

24.6 

2.6 

27.2 

5.0 km 

9.6 

1.05 

10.65 

7. 5 km 

5.8 

0.34 

6.14 

10. 0 km 

3.53 

0.12 

3.65 

(b) AFCRL Midlatitude Winter Clear Atmosphere and C^ = 

4.2 X 10”^^ 


ilV 


17.3 

3.4 

20.7 


10.5 

0. 3 

10.8 


6 . 5 

- 

6 . 5 

7 km 

3.5 

- 

3. 5 

10 Ion 

1.8 


1.8 


VII. CONCLUSIONS 


The attenuation of the CO 2 laser radiation due to absorption, scattering, 
and turbulence of the atmosphere is estimated for the January 1973 flight test 
conditions. The theoretical values of S/N including the atmospheric attenua- 
tion compare favorably with the flight test results. 
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The upper and lower limits on the atmospheric attenuation at each 
altitude up to 10 km may be read from Figure 6 for a two-way horizontal path 
of 20 km range. These calculations are based on the currently accepted models 
of atmospheric pressure, temperature, humidity, aerosol concentration, and 
turbulence strength. These results are expected to be useful in the interpreta- 
tion of the future test results at several altitudes using the pulsed CO 2 laser 
Doppler system. 
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